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Abstract 

We consider a D7-brane probe of AdSs x S 5 in the presence of pure gauge B- 
field. The dual gauge theory is flavored Yang-Mills theory in external magnetic 
field. We explore the dependence of the fermionic condensate on the bare quark 
mass m q and study the discrete self-similar behavior of the theory near the 
origin of the parametric space. We calculate the critical exponents of the bare 
quark mass and the fermionic condensate. A study of the meson spectrum 
supports the expectation based on thermodynamic considerations that at zero 
bare quark mass the stable phase of the theory is a chiral symmetry breaking 
one. Our study reveals the self-similar structure of the spectrum near the 
critical phase of the theory, characterized by zero fermionic condensate and we 
calculate the corresponding critical exponent of the meson spectrum. 



1 Introduction 



The existence of gauge/sting dualities has been anticipated from the very beginning 
of these fields, based on the striking similarity between the large N t'Hooft limit 
and the genus expansion in string theory. Further insights were revealed with the 
discovery of D-branes, their description in terms of DBI action and identification as 
sources of the known p-brane solutions in super gravity. However it was not until 
Maldacena conjectured his decoupling limit when significant progress was made and 
the AdS/CFT correspondence emerged p] as a manifestation of the gauge/string 
duality providing holographic description of super gravity on AdS 5 x S 5 space in 
terms of J\f—4 SUSY Yang-Mills theory living at the asymptotic boundary. The 
observed correspondence was conjectured to hold for the full string theory on any 
asymptotically AdS§ x S 5 background. One of the most remarkable features of the 
correspondence is that it is a strong-weak correspondence and that it can give us 
tools to explore the strongly coupled regimes of the Yang-Mills theory. 

In recent years progress has been made towards the study of matter in fundamen- 
tal representation in the context of AdS/CFT correspondence. One way to achieve 
this is by introducing space filling flavor D7-branes in the probe limit [2] and in order 
to keep the probe limit valid the condition Nf <C N c is imposed. The fundamental 
strings stretched between the stack of N c D3 branes and the flavor Nf D7-branes give 
rise to J\f=2 hypermultiplet, the separation of the D3 and D7 branes in the transverse 
directions corresponds to the mass of the hypermultiplet, the classical shape of the 
D7-brane encodes the value of the fermionic condensate and its quantum fluctuations 
describe the light meson spectrum of the theory [3]. This technique for introducing 
fundamental matter has been widely employed in different backgrounds. Of particu- 
lar interest was the study of non supersymmetric backgrounds and phenomena such 
as spontaneous chiral symmetry breaking. These phenomena were first studied in 
this context in pE], where the authors developed an appropriate numerical technique. 
In recent years this approach received further development, and has proven itself as 
powerful tool for the exploration of confining gauge theories, in particular, for the de- 
scription of their thermodynamic properties or for the building of phenomenological 
models relevant to Q CD 0- [IS]. 

The paper is organized as follows: 
In the second section we review the method of introducing magnetic field to the 
theory, employed in [9]. We describe the basic properties of the D7 brane embed- 
ding and the thermodynamic properties of the dual gauge theory, in particular the 
dependence of the fermionic condensate on the bare quark mass. We describe the 
spontaneous chiral symmetry breaking caused by the external magnetic field and 
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comment on the spiral structure in the condensate vs. bare quark mass diagram. 

The third section contains our main results and splits into two parts: 
The first part is dedicated to the detailed study of the spiral structure described 
in [9]. We perform analysis similar to the one considered in [5] for the study of 
merger transitions and calculate the critical exponents of the bare quark mass and 
the fermionic condensate. We also describe the discrete self-similarity of the spiral 
and calculate the scaling factor characterizing it. 

In the second part of this section we consider the meson spectrum of the states 
corresponding to the spiral. First we study the critical embedding corresponding 
to the center of the spiral and reveal an infinite tower of tachyonic states organized 
in a decreasing geometrical series. Next we consider the dependence of the meson 
spectrum on the bare quark mass and confirm the expectations based on thermody- 
namic considerations that only the lowest branch of the spiral is stable. We observe 
that at each turn of the spiral there is one new tachyonic state. We comment on 
the self-similar structure of the spectrum and calculate the critical exponent of the 
meson mass. We also consider the spectrum corresponding to the lowest branch of 
the spiral and for a large bare quark mass reproduce the result for pure M = 2 
Supersymmetric Yang Mills Theory obtained in [3]. 

We end with a short discussion of our results and the possible directions of future 
study. 

2 Fundamental matter in external magnetic field 

In this section we briefly review the method of introducing external magnetic field 
to the theory considered in [9] and the basic properties of the D7-brane probe in this 
background. We also review the properties of the corresponding dual theory and the 
effect that the external magnetic field has on it. 
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2.1 Basic Configuration 

Let us consider the AdS?, x S 5 geometry describing the near horizon geometry of a 
stack of N c extremal D3-branes. 

ds 2 = t-A-K + dx 2 ) + R 2 —r + R 2 dttl (1) 



R 2 " u ' v? 

4 



9sC(4) = J^dx A dx 1 A dx 2 A dx 3 , 
R 4 = ^g s N c a' 2 . 

In order to introduce fundamental matter we first rewrite the metric in the fol- 
lowing form : 

2 _|_ jJ2 

ds 2 = — [-dx 2 + dx\ + dx\ + dx 2 ] + ^ - ^ 2 [dp 2 + p 2 dQ 2 3 + dL 2 + L 2 d<p 2 ] , 

dtt 2 3 = dip 2 + cos 2 ipd(3 2 + sin 2 ipd-f 2 , (2) 

where p, if), (5, 7 and L, are polar coordinates in the transverse 1Z 4 and 7£ 2 respec- 
tively, satisfying: u 2 = p 2 +L 2 . Next we use £0,1,2,3, Pi V*j A 7 to parametrise the world 
volume of the D7-brane and consider the following ansatz [2] for it's embedding: 

= const, (3) 
L = L(p) . 

Leading to the following form of the induced metric: 

ds = p2 + ^ (P)2 [-dx 2 + tfc 2 + dx\ + tfc 2 ] + p2 * [(1 + L'(p) W + p 2 dQl) . (4) 
Now let us consider the NS part of the general DBI action: 

Sdbi = -— [ d^det l ' 2 {P[G ab + B ab ] + 2na'F ab ) . (5) 

Ms 

Here pi = [(27r) 7 )a /4 ] _1 is the D7-brane tension, P[G ab ] and P[B ab ] are the in- 
duced metric and -B-field on the D7-brane's world volume, while F ab is its gauge field. 
A simple way to introduce magnetic field would be to consider pure gauge B-field 
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along the "flat" directions of the geometry xq — x% corresponding to the D3-branes 
world volume: 

B {2) = Hdx 2 A dx 3 . (6) 

The constant H is proportional to the magnetic component of the EM field. Note 
that since the 5-field is a pure gauge dB = the corresponding background is still a 
solution to the supergravity equations. On the other side the gauge field F ab comes 
in next order in a' expansion compared to the metric and the -B-field components. 
Therefore to study the classical embedding of the D-brane one can leave only the 
G a b + B ao part of the DBI-action. It was argued in [9] that one can consistently 
satisfy the constraints imposed on the classical embedding resulting from integrating 
out F a b- The resulting effective lagrangian is: 



£ = -^p 3 sin^cosVVl + L' 2 Jl + — — . (7) 

9s V {p 2 + L 1 ) 2 

The equation of motion for the profile Lo(p) of the D7-brane is given by: 

2p 3 L R 4 H 2 




mi (P 2 + Llf 



. 



As expected for large (Lq + p 2 ) — > oo or H —> 0, we get the equation for the pure 
AdS 5 x S 5 background [2]: 

8 P ( p 3 , L '° 1 = . 
Therefore the solutions to equation (J8T) have the following behavior at infinity: 



L Q 0) = m + 4 + • • • , (9) 
P 2 

where the parameters m (the asymptotic separation of the D7- and D3- branes) 
and c (the degree of bending of the D7-brane) are related to the bare quark mass 
m q = m/2na' and the fermionic condensate (ipip) oc — c respectively [6]. We have 
provided derivation of these relations in Appendix A. As we shall see below, the 
presence of the external magnetic field and its effect on the dual SYM provide a non 
vanishing value for the fermionic condensate, furthermore the theory exhibits chiral 
symmetry breaking. 
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Now notice that H enters in (J7| only through the combination H 2 R 4 . The other 
natural scale is the asymptotic separation m. It turns out that different physical 
configurations can be studied in terms of the ratio m 2 = m 2 /(HR 2 ): Once the m 
dependence of our solutions are known, the m and H dependence follows. Indeed let 
us introduce dimensionless variables via: 

p = RVHp , L = R^/HL , L' {p) = L\p) . (10) 

The equation of motion d8j) then takes the form: 



d- P p 3 , — = J i + — - + i — = o . ii 

1 + V (P 2 + L 2 ) 2 ) Ji + w A z _(p 2 + L 2 r 
The solutions for L(p) can be expanded again to: 

L(p)=m + ^ + ..., (12) 
P 2 



and using the transformation (10) we can get: 

c = dR 3 H 3/2 . (13) 

2.2 Properties of the Solution 

The properties of the solution have been explored in [9], both numerically and ana- 
lytically, when possible. Let us briefly review the main results. 

For weak magnetic field H and non-zero bare quark mass m it was shown that 
the theory develops a fermionic condensate: 

mcx-c = -^H 2 , (14) 
or using dimensionless variables: 

Am 

The case of strong magnetic field H can be explored by numerically solving 



equation (11), it is convenient to use initial conditions in the IR as has been recently 
discussed in the literature [7], [8]. We used the boundary condition L'(p)\p =0 = 0. 
We used shooting techniques to generate the embedding of the D7 for a wide range 
of m. Having done so we expanded numerically the solutions for L(p) as in equation 
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(12) and generated the points in the (m, — c) plane corresponding to the solutions. 
The resulting plot is presented in figure [T} 

As one can see there is a non zero fermionic condensate for zero bare quark mass 
and hence there is a Spontaneous Breaking of the Chiral Symmetry. The corre- 
sponding value of the condensate is c cr = 0.226. It is also evident that the analytical 
expression for the condensate (15) that we got in the previous section is valid for 
large in, as expected. Now using equation (13) we can deduce the dependence of c cr 
on H: 

c cr = c cr R 3 H 3/2 = 0.22QR 3 H 3/2 . (16) 

Another interesting feature of our phase diagram is the spiral behavior near the 
origin of the (m, — c)-plane which can be seen in figure [3] Note that the spiral 
presented in this figure has two arms, we have used the fact that any two points 
in the (m, — c) plane related by reflection with respect to the origin describe the 
same physical state. A similar spiraling feature has been observed in ref. [8], where 
the authors have argued that only the lowest branch of the spiral corresponding to 
positive values of m is the stable one (corresponding to the lowest energy state). The 
spiral behavior near the origin signals instability of the embedding corresponding to 
Lq = 0. If we trace the curve of the diagram in figure [3] starting from large m, 
as we go to smaller values of m we will reach zero bare quark mass for some large 
negative value of the fermionic condensate c cr . Now if we continue tracing along the 
diagram one can verify numerically that all other points correspond to embeddings 
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of the D7-brane which intersect the origin of the transverse plane at least once. 
After further study of the right arm of the spiral, one finds that the part of the 
diagram corresponding to negative values of in represents solutions for the D7-brane 
embedding which intersect the origin of the transverse plane odd number of times, 
while the positive part of the spiral represents solutions which intersect the origin of 
the transverse plane even number of times. The lowest positive branch corresponds 
to solutions which don't intersect the origin of the transverse plane and is the stable 
one, while the upper branches have correspondingly 2,4, etc., intersection points and 
are ruled out after evaluation of the free energy. Indeed let us explore the stability of 
the spiral by calculating the regularized free energy of the system. We identify the 
free energy of the dual gauge theory |3B], [SU] with the wick rotated and regularized 
on-shell action of the D7-brane: 

F = 27r 2 N f T D7 R i H 2 I D7 , (17) 

pmax / I \ 

I m = j dp [p^l + — | , a) Vl + L* - 1 (18) 



The second term under the sign of the integral in (18), corresponds to the subtracted 



free energy of the L(p) = embedding and serves as a regulator. Now we can 



evaluate numerically the integral in (18) for the first several branches of the spiral. 
The corresponding plot is presented in figure ^ Note that we have plotted Id7 versus 
|m|, since the bare quark mass depends only on the absolute value of the parameter 
rh. The lowest curve on the plot corresponds to the lowest positive branch of the 
spiral, as one can see it has the lowest energy and thus corresponds to the stable 
phase of the theory. 

In the next section we will provide more detailed analysis of the spiral structure 
from Figure [3] and explore the discrete self-similarity associated to it. 

3 Criticality and Spontaneous chiral symmetry break- 
ing 

3.1 The Spiral Revisited 

In the following section we analyze the spiral structure described in [Q]. The technique 
that we employ is similar to the one used in [5] and [28] , where the authors studied 
merger transitions in brane-black-hole systems. 
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Figure 2: The lowest lying curve correspond to the positive rh part of the lowest branch 
of the spiral, suggesting that this is the stable phase of the theory. 



Let us explore the asymptotic form of the equation of motion of the D7-brane 
probe (11 ) in the near horizon limit p 2 + L 2 — > 0. To this end we change coordinates 
to: 



p — ► Ap; L — > XL; 
and consider the limit A — > 0. The resulting equation of motion is: 

V 



P 



^P 2 + L 2 Vl + L' 2 



+ 2\ 1 + L' 2 - 



p L 



(p 2 + L 2 ) 2 







Equation (20) enjoys the scaling symmetry: 

P -»• m L-> /iL; 



(19) 



(20) 



(21) 



In the sense that if L = f(p) is a solution to the E.O.M. then -f(p:p) is also a 
solution. Next we focus on the region of the parametric space, close to the trivial 
L = embedding, by considering the expansion: 



L = + (2vra')x 
and linearizing the E.O.M. . The resulting equation of motion is: 

Pdp(pdpx) + 2x = 

and has the solution : 

X = Acos(v / 21np) + J Bsin(\ / 21np) . 



(22) 
(23) 
(24) 
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Figure 3: A magnification of figure [T] to show the spiral behavior near the origin of 
the (— c, m)-plane. We have added the second (left) arm of the spiral representing the 
(m, — c) — ► (—77t, c) symmetry of the diagram. 



(25) 



Now under the scaling symmetry p — ► pp the constants of integration A and B 
transform as: _ 

A \ 1 / cos v^ln /i sin v^ln/A (A" 

B J ^ \ — sin \/21n /i cos-\/21n/i/ \-B 

The above transformaton defines a class of solutions represented by a logarithmic 
spiral in the parametric space (A, B) generated by some (A in , B in ), the fact that we 
have a discrete symmetry x — X suggests that (—A in , —B in ) is also a solution and 
therefore the curve of solutions in the parametric space is a double spiral symmetric 
with respect to the origin. Actually as we are going to show there is a linear map 
from the parametric space (A, B) to the plane (m, — c) which explains the spiral 
structure, a subject of our study. To show this let us consider the linearized E.O.M. 
before taking the A — > limit : 

py/l+ffidpipy/l + pdtx) + 2x = 0; x = Ax; , (26) 
with the solution: 

icosy^ln - 9 = + Bsiny^ln , 9 . (27) 



y 

'i + \A + p 4 Ji + ^i + p 



7A 



Expanding at infinity: 



c t B 1 

X = m+ — + --- = A-— — + ..., 28 

P 2 V2 P 
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we get: 



■bU) ■ (29) 



Now if we match our solution (27) with the solution in the p — > limit (24) we 
should identify (A, B) with the parameters (A, B). Combining the rescaling property 
of (A, B) with the linear map to (m, — c) we get that the embeddings close to the 
trivial embedding L = are represented in the (m, — c) plane by a double spiral 
defined via the transformation: 



m 




a/2 In p — \f2 sin \/2\np 
sin y/2 In p cos y/2 In p 




(30) 



Note that the spiral is double, because we have the symmetry (m, — c) — > (— rh,c). 
This implies that in order to have similar configurations at scales /ii and pi we should 
have: 

cos \pl In /i x = ± cos v^ln/xs (31) 

and hence : 

v^ln^^-mr, (32) 
Pi 

which is equivalent to: 

Therefore we obtain that the discrete self-similarity is described by a rescaling by a 
factor of: 

q = e' w/ ^ ss 0.10845 . (34) 

This number will appear in the next subsection where we will study the meson 
spectrum. As one may expect the meson spectrum also has a self-similar structure. 

It is interesting to confirm numerically the self-similar structure of the spiral 
and to calculate the critical exponents of the bare quark mass and the fermionic 
condensate. It is convenient to use the separation of the D3 and D7 branes at p = 0, 
L in = 1/(0) as an order parameter. There is a discrete set of initial separations L in , 
corresponding to the points H , Hi, H 2 , . . . in figure [3] , for which the corresponding 
D7 brane's embeddings asymptote to m — = as p —>■ oo. The trivial L = 
embedding has Lj n = and is the only one which has a zero fermionic condensate 
(c = 0), the rest of the states have a non zero c and hence a chiral symmetry is 
spontaneously broken. Each such point determines separate branch of the spiral 
where c = c(m) is a single valued function. On the other side each such branch 
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has both positive fh and negative fh parts. The symmetry of the double spiral from 
figure [3j suggests that the states with negative fh are equivalent to positive fh states 
but with an opposite sign of c. This implies that the positive and negative fh parts of 
each branch correspond to two different phases of the theory, with opposite signs of 
the condensate. As we can see from figure [2] the lowest positive branch of the spiral 
has the lowest free energy and thus corresponds to the stable phase of the theory. 
In the next subsection we will analyze the stability of the spiral further by studying 
the light meson spectrum of the theory near the critical L = embedding. 

Here we are going to show that both the bare quark mass fh and the fermionic 
condensate c have critical exponent one, as L in — > 0. Indeed let us consider the scaling 
property (25), (30). If we start from some V\ n and transform to L in = -L® n , we can 
solve for /i and using equation (30) we can verify that the bare quark mass and the 
fermionic condensate approach zero linearly as Lj n — > 0. To verify numerically our 
analysis we generated plots of fh/L in vs. \[2 logZj n /27r and c/L in vs. \[2 \og L in /2n 
presented in figure |4j 

The red curves in these figures represent a fit with trigonometric functions of a 
unit period, as one can see the fit is very good as L in — > 0. On the other side for 
large L in we obtain the results for a pure AdSs x S 5 space, namely L = const, c = 0. 
It is also evident from the plots that the critical exponents of fh and c are equal to 
one. 



3.2 The Meson Spectrum 

In this section we will explore the light meson spectrum of the theory corresponding 
to quadratic fluctuations of the D7 brane embedding. In particular we will consider 
the spectrum corresponding to the fluctuations of L. The equations of motion of 
the fluctuation modes were derived in [U] and it was shown that the vector and the 
scalar spectrum mix due to the non-zero magnetic field. Some interesting effects 
such as Zeeman splitting of the states and a characteristic y/m dependence of the 
meson spectrum have been reported. However the analysis performed in [9] is only 
for the fluctuations along 0, for the lowest positive branch of the spiral from figure 
[3] (the one corresponding to point H ). In this letter we extend the analysis of 
the spectrum to all branches of the spiral (points H 1 , H 2 , . . . in figure [3]) and show 
that the ground states of all inner branches of the spiral are tachyonic, proving 
that the phases described by these branches of the spiral are unstable as opposed 
to metastable. Our analysis reveals the self-similar structure of the spectrum and 
we obtain the critical exponents of the tachyonic spectrum as one approaches the 
critical L = embedding. The chapter is organized as follows: 
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1.01- 




-0.6 



Figure 4: The red curves represent fit with trigonometric functions of unit period. For 
small Li n the fit is very good, while for large Zj n we get the results for pure AdS^ x S 5 
space, namely L = const, c = 0. The plots also verify that the critical exponents of rh and 
c are equal to one. 
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First we study the spectrum of the L = embedding in the spirit of the anal- 
ysis provided in [32] ■ We perform both a numerical and analytical study and show 
that the spectrum contains infinitely many tachyonic states approaching zero in a 
decreasing geometrical series, representing the self-similar structure of the meson 
spectrum. 

Next we study the spectrum as a function of the bare quark mass and show that 
at each turn of the spiral one of the energy levels become tachyonic. Similar behavior 
has been recently reported in [33] • We show that as we approach the critical L = 
embedding the spectrum becomes tachyonic and the corresponding critical exponent 
is two. We also present plots showing the spiraling of the spectrum as one approaches 
criticality. 

Finally we provide an analysis of the spectrum of the stable branch of the spiral 
and comment on the small m behavior of the spectrum as a consistent with the 
spontaneous chiral symmetry breaking scenario. 

3.2.1 The critical L = embedding 

In this section we study the L = embedding and in particular the spectrum of the 
fluctuations along the L coordinate. Let us go back to dimensionfull coordinates and 
consider the following change of coordinates in the transverse i? 6 space: 



In these coordinates the trivial embedding corresponds to 9 = and in order to 
study the quadratic fluctuations we perform the expansion: 



Note that in order to study the mass spectrum we restrict the D7 brane to fluc- 
tuate only in time. In a sense this corresponds to going to the rest frame. Note that 
due to the presence of the magnetic field there is a coupling of the scalar spectrum to 
the vector one, however for the fluctuations along 9 the coupling depends on the mo- 
menta in the (x2, X3) plane and this is why considering the rest frame is particularly 
convenient . 

Our analysis follows closely the one considered in [32], where the authors have 
calculated the quasinormal modes of the D7-brane embedding in the AdS-black hole 



p = u cos 9 , 
L = u sin 9 . 



(35) 




(36) 
(37) 



56 = e~ int h(u) . 



13 



background by imposing an in-going boundary condition at the horizon of the black 
hole. Our case is the T — > limit and the horizon is extremal, however the 6 = 
embedding can still have quasinormal excitations with imaginary frequencies, corre- 
sponding to a real wave function so that there is no flux of particles falling into the 
zero temperature horizon. The resulting equation of motion is: 

h+ {u + ^Tmp) h + {^ + u*) h = - (38) 

It is convenient to introduce the following dimensionless quantities: 

R % - = ^; , (39) 



and make the substitution jl2] 



leading to the equation for the new variable f(z): 

f\z) + {u? - V(z)) f(z) = . (41) 
Where the effective potential is equal to: 



The potential in (42) goes as ^ for z —* and as — for z — >• oo and is presented 
in figure [5j As it was discussed in |32] if the potential gets negative the imaginary 
part of the frequency may become negative. Furthermore the shape of the potential 
suggests that there might be bound states with a negative us 2 . To obtain the spectrum 



we look for regular solutions of (41) imposing an in-falling boundary condition at 
the horizon (z — > oo). 

The asymptotic form of the equation of motion at z — > oo is that of the harmonic 
oscillator: 

f"(z)+u J 2 f(z) = 0, (43) 

with the solutions e ±luJZ , the in-falling boundary condition implies that we should 
choose the positive sign. In our case the corresponding spectrum turns out to be 
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V(z) 




Figure 5: A plot of the effective potential V(z) given in equation (42) 



tachyonic and hence the exponents are real. Therefore the in-falling boundary condi- 
tion simply means that we have selected the regular solution at the horizon: z — > oo. 
We look for a solution of the form: 



f(z) = e + ^ z S(z) . 
The resulting equation of motion for S(z) is: 

(-3 - Qz 4 + 9z 8 )S(z) + 4z 2 (l + z 4 f (2iuS'(z) + S"(z)) = . 



(44) 



(45) 



Next we study numerically equation (45). After solving the asymptotic form of 
the equation at the Horizon, we impose the following boundary condition at z = 1/e, 
where e is a numerically small number typically e = 10~ 9 : 



S(l/e) = 1 - 



9ie 



S\l/e) 



9ie^ 
8a; 



(46) 



after that we explore the solution for a wide range of u — iu>i. We look for regular 
solutions which have |5(e)| ~ 0, this condition follows from the requirement that 
X oc z 3 as z — > 0. It turns out that regular solutions exist for a discrete set of 
positive w;<l. The result for the first six modes that we obtained is presented in 
table CD 

The data suggests that as ujj — > the states organize in a decreasing geometrical 
series with a factor q w 0.1084. Up to four significant digits, this is the number 



from equation (34), which determines the period of the spiral. We can show this 
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Table 1: 



n 




. An) / (n-1) 





2.6448 x 10" 1 




1 


2.8902 x 1(T 2 


0.10928 


2 


3.1348 x 10~ 3 


0.10846 


3 


3.3995 x 1(T 4 


0.10845 


4 


3.6865 x 1(T 5 


0.10844 


5 


3.9967 x 10~ 6 


0.10841 



analytically. To this end let us consider the rescaling of the variables in equation 
(45) given by: 

z = Xz; lj = u/\; X — > oo; . (47) 

This is leading to: 

9S{z) + Az 2 (2iCjS\z) + S"{z)) + 0(\- 4 ) = . (48) 



The solution consistent with the initial conditions at infinity (46) can be found to 
be: 

S{z) = ] ^e- l ^e- i ^V^H (1 Uzuy, Co = iuf, , (49) 

2 i\l 2, 

where is the Hankel function of the first kind. Our next assumption is that in 
the uji — > limit, this asymptotic form of the equation describes well enough the 
spectrum. To quantize the spectrum we consider some Zq — Zq/X <C 1, where we have 
1 <C zq <C A so that the simplified form of equation (48) is applicable and impose: 

3(zq) = . (50) 

Using that zCj = izuji this boils down to: 

H^iujjzo) = . (51) 

Now using that ujiZq <C 1 for a sufficiently small Ui, we can make the expansion: 

H^(icujz ) « -At (W )^ - (ujz )- lV ~ 2 ) +iA 2 (Wo)^ + M- iV§ ) , (52) 
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where Ai and A 2 are real numbers defined via: 



A x + iA 2 



—i(i/2)- 

71 



iV2 



T(iV2) . 



(53) 



This boils down to: 



cos 



(\/21n( 



UiZq) + 



The first equation in (54) leads to: 



(n) 



0: 



V2 e ~ 



= n/2~aig(A 1 + iA 2 ); 



'—2 



z 



(0) 

uj) 'q 



suggesting that: 



q = e ^ » 0.10845 



(54) 



(55) 



(56) 



This is the number given in (34). Note that the value of zn is a free parameter that 
we can fix by matching equation (55) to the data in table [lj On the other side S(z) 
given in equation (49) depends only on zuj = icojz and therefore once we have fixed 



zq we are left with a function of cu/, which zeroes determine the spectrum, equation 
(50). It is interesting to compare it to the numerically obtained plot of \S(e)\ vs. Ui, 



that we have used to determine the spectrum numerically. The result is presented in 
figure [6j where we have used the n = 3 entry from table [T] to fix Zq. One can see the 
good agreement between the spectrum determined by equation (50), the red curve 



in figure pi and the numerically determined one, the dotted blue curve. 



3.2.2 The Spectrum near criticality 

In this chapter we study the light meson spectrum of the states forming the spiral 
structure in the (m, — c) plane, figure [3j In particular we focus on the study of 
the fluctuations along L. The corresponding equation of motion was derived in 
The effect of the magnetic field H is to mix the vector and the meson parts of the 
spectrum. However if we consider the rest frame by allowing the fluctuations to 
depend only on the time direction of the D3 branes' world volume, the equation of 
motion for the fluctuations along L decouple from the vector spectrum. To this end 
we expand: 

L = L (p) + (27ra') X (p,t) , (57) 
X = h(p) cos Mt . 
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|S(e)| 



|S(e)| vs log(w/) 




-3.0 



-2.5 



-2.0 



-1.5 



-1.0 



-0.5 



log [to/] 



Figure 6: The dotted blue curve corresponds to the numerical solution of equation (45), 



while the thick red curve is the one determined by equation (50). The plots are scaled to 
match along the vertical axis. 



Here L (p) is the profile of the D7 brane's classical embedding. The resulting equa- 
tion of motion for h(p) is: 



d 2 9 



dg L 'o 



) 2 1+L' 2 dL 2 ^ W P\8L l+Ltf 



h = 



where 



,(p,L ,L' )=p 3 v /l + L '\ 1 + 



5V) ■ ' 'I' V ■' ' -M! \i ' T ( p 2 +L 2)2 • 

It is convenient to introduce the dimensionless variables: 

P 



h 



h 



RVH 



f L „ 

j -^0 — 7=] P 



RVH 



M 



MR 



H 



leading to: 



with 



hi 



+ 9 



M 2 



_ £i + q. 



L'n 



(p2 + qyi+L' 2 8L 2 1 "P\8L l+L% 



g(p,L ,L' ) = p\/l + L' n 2 Jl + 



(58) 



h = 



(59) 



(60) 



(P 2 +Ll) 2 • 

We study the normal modes of the D7 brane described by equation (60) by 



imposing Neumann boundary conditions at p = 0. Since our analysis is numerical 
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we solve the equation of motion (60) in terms of a power series for small p and 
impose the appropriate initial conditions for the numerical solution at p — e, where e 
is some very small number. In order to quantize the spectrum we look for numerical 
solutions which are normalizable and go as 1/ p 2 at infinity. 

Let us study the dependence of the spectrum of M on the bare quark mass rh, 
for the states corresponding to the spiral structure from figure [3] A plot of the 
spectrum of the first three excited states is presented in figure [7j The classification 
of the states in terms of the quantum number n is justified, because at large rh the 
equation of motion for the fluctuations asymptotes to the equation of motion for the 
pure AdS 5 x S 5 space, considered in [3], where the authors obtained the spectrum in 
a closed form. Note that the diagram has a left-right symmetry. This is because we 
plotted the spectrum for both arms of the spiral in order to emphasize its self-similar 
structure, physically only one side of the diagram is sufficient. 

Let us trace the blue curve corresponding to the n = state starting from the 
right-hand side. As rh decreases the mass of the meson decreases and at rh — it 
has some non-zero value. This part of the diagram corresponds to the lowest positive 
branch of the spiral from figure [3] (the vicinity of point H ). It is satisfying to see 
that the lowest positive rh branch of the spiral is tachyon free and therefore stable 
under quantum fluctuations. Note that despite that the negative rh part of the lowest 
branch has no tachyonic modes in its fluctuations along L, it has a higher free energy 
(as can be seen from figure [2]) and is thus at best metastable. 

One can also see that the spectrum drops to a zero and becomes tachyonic exactly 
at the point where we start exploring the upper branch of the spiral. This proves 
that all inner branches correspond to true instability of the theory and cannot be 
reached by super-cooling. As we go deeper into the spiral, the n = spectrum 
remains tachyonic and spirals to some critical value. The dashed line denoted by 
uf^ in figure [i] corresponds to the first entry in table 1 As one can see this is the 
critical value approached by the spectrum. 

Now let us comment on the n = 1, 2 levels of the spectrum represented by the red 
and green curves, respectively. As one can see the n = 1 spectrum becomes tachyonic 
when we reach the third branch of the spiral (the vicinity of point H 2 in figure [3]) 
and after that follows the same pattern as the n = level, spiraling to the second 
entry oo^ from table 1 The n = 2 level has a similar behavior, but it becomes 
tachyonic at the next turn of the spiral and it approaches the next entry from table 
[T] Similar feature was reported recently in jH] where the authors studied topology 
changing transitions. The above analysis suggests that at each turn of the spiral, 
there is one new tachyonic state appearing. It also suggests that the structure of the 
n-th level is similar to the structure of the n + 1-th level and in the n —>■ oo limit 
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-0.10 



-0.15 




-0.30 



Figure 7: A plot of the meson spectrum corresponding to the two arms of the spiral 
structure at the origin of the (m, — c) plane. The ground state (n = 0) becomes tachyonic 
for the inner branches of the spiral, while only the lowest branch is a tachyon free one. The 
tachyon sector of the diagram reveals the self-similar structure of the spectrum. 



20 



this similarity becomes an exact discrete self- similarity. The last feature is apparent 
from the tachyonic sector of the diagram in the second plot in figure [7j the blue, red 
and green curves are related by an approximate scaling symmetry, the analysis of 
the spectrum of the critical L = embedding suggests that this symmetry becomes 



exact in the n — > oo limit with a scaling factor of q given in equation (34). 

It is interesting to analyze the way the meson mass M approaches its critical value 
and compute the corresponding critical exponent. Let us denote the critical value of 
M by M* and consider the bare quark mass rh as an order parameter, denoting its 
critical value by m*. We are interested in calculating the critical exponent a defined 
by: 

\M-M*\ oc \rh-rh*\ a . (61) 

We will provide a somewhat heuristic argument that a = 2 and will confirm this 
numerically. To begin with let us consider the energy density of the gauge theory 
E as a function of the bare quark mass rh. Now let us consider a state close to the 
critical one, characterized by: 

M .(/. • 5M] rh = m* + 5rh; E = E* + 5E; . (62) 

Next we assume that as we approach criticality the variation of E and M are pro- 
portional to the variation of the energy scale and hence 8E oc 8M. Therefore we 
have: 

^Ac, (63) 
dm dm 



where c is the fermionic condensate. The second relation in (63) was argued in [6]. 
In the previous section we argued that the critical exponent of the condensate is one 
and since the critical embedding has a zero condensate it follows that c oc \m — m*|. 
Therefore we have: 

5M . . „ „ . , , 

— - oc a \m — mJ oc \m — mJ (64) 
dm 

and hence a = 2. 

Now let us go back to figure [7j As we discussed above, for each energy level n 
the tachyonic spectrum spirals to the critical value ui^, corresponding to the center 
of the spiral. If we focus on the rh = axis, we can see that for each level we have a 
tower of tachyonic states at a zero bare quark mass, corresponding to the different 
branches of the spiral. Let us denote by the imaginary part of the meson 

spectrum, corresponding to the fc-th tachyonic state of the n-th energy level, at a 
zero bare quark mass rh. As we go deeper into the spiral, k — ► oo and — > Mi n \ 
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the data in figure 7 suggests that M, 



r(n) 



UJ 



(n) 



On the other side if the meson 



spectrum has a critical exponent of two, one can show that for a large k: 



M, 



(n) 



M 



(n) 



M 



(n) 
fc-1 



Ml 



(») 



(65) 



where q is given by equation (34). We can solve for M* 



Mi n) = M, 



fc-i 



+ 



1-g 2 



(66) 



Now assuming that for k 



1,2 the approximate geometrical series defined via 
(65) is already exact we calculate numerically Af}"^ , M^ for the n = 0,1,2 levels 
and compare the value of M^ obtained by equation (66) to the first three entries 
in table [TJ The results are presented in table [2} 



Table 2: 



n 


M[ n) 


M^ n) 


Mi n) 


.in) 





2.7530 x HT 1 


2.6460 x 10" 1 


2.6447 x 10" 1 


2.6448 x 10" 1 


1 


3.0162 x 10" 2 


2.8917 x 10~ 2 


2.8902 x 10" 2 


2.8902 x 10" 2 


2 


3.2715 x 10" 3 


3.1363 x 10" 3 


3.1347 x 10" 3 


3.1348 x 10" 3 



One can see that up to four significant digits the critical value of the meson 
spectrum is given by the imaginary part of the quasi normal modes presented in 
table [T] This supports the above argument that the meson spectrum has a critical 
exponent of two. Another way to justify this, is to generate a plot of the meson 
spectrum similar to the one presented in figure [4] for the bare quark mass m and the 
fermionic spectrum c. Notice that M approaches criticality from above, while the 
parameter m oscillates around the critical value = 0. This suggests to use M as 
an order parameter and to generate a plot of rh/(M — M*) 2 vs. v^21og \ M — M*\/2k. 



Note that according to equation (|65|) the plot should represent periodic function of 

as a 



an unit period. The resulting plot for the n = level, using Ml ** from table 



critical value, is presented in figure M 



3.2.3 The stable branch of the spiral 



In this subsection we consider the spectrum corresponding to the states far from the 
origin of the (m, — c), which is the outermost branch of the spiral ending at point 
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m 




Figure 8: A plot of the bare quark mass meson vs. the meson spectrum, in an appropriate 
parameterization, determined by the critical exponents of fh and M. The discrete self- 
similar structure of the spectrum is manifested by the periodicity of the plotted function. 



Hq from figure [3j The fluctuations of the D7-brane corresponding to the massless 
scalar were studied in [9] and some features consistent with the spontaneous chiral 
symmetry breaking, such as a characteristic yfm behavior [51] were reported. 

Here we complement the analysis by presenting the results for the fluctuations 
along the L coordinate. Since this is the massive field in the spontaneous chiral sym- 
metry breaking scenario, we expect a \J const + fh behavior of the meson spectrum 
for small values of fh. Note that such a behavior simply means that the spectrum of 
the L fluctuations has a mass gap at zero bare quark mass and that the slope of the 
spectrum vs. the bare quark mass function is finite. It is satisfying that our results 
are in accord with this expectations. 



To obtain the spectrum, we solve numerically equation (60) imposing Neumann 
boundary conditions at p = 0. A plot of the first five energy levels is presented 
in figure |9j As one can see at large fh the spectrum approximates that of the pure 
M = 2 Flavored Yang Mills theory studied in [3J, where the dependence of the meson 
spectrum on the bare quark mass was obtained in a closed form: 

2m 



M = — y/( n + l + l)(n + l + 2) . (67) 

Here / is the quantum number corresponding to the angular modes along the internal 
S 3 sphere wrapped by the D7 brane and is zero in our case. After introducing the 
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dimensionless variables defined in (59), equation (67) boils down to: 

M = 2y/(n + l)(n + 2)m . 



(68) 



The black dashed lines in figure 9 represent equation ( 68 ) . The fact that the meson 



spectrum asymptotes to the one described by (68) justifies the use of the quantum 
number n to classify the meson spectrum. One can also see that as expected the 
spectrum at zero bare quark mass has a mass gap. 




Figure 9: A plot of the meson spectrum corresponding to the stable branch of the spiral. 



The black dashed lines correspond to equation ( 68 ) , one can see that for large rh the meson 
spectrum asymptotes to the result for pure AdSs x S 5 space. One can also see that at zero 
bare quark mass rh there is a mass gap in the spectrum. 



4 Conclusion 

In this paper we performed a detailed analysis of the spiral structure at the origin of 
the condensate vs. bare quark mass diagram. We revealed the discrete self-similar 
behavior of the theory near criticality and calculated the corresponding critical ex- 
ponents for the bare quark mass, the fermionic condensate and the meson spectrum. 

Our study of the meson spectrum confirmed the expectations based on thermo- 
dynamic considerations that the lowest positive rh branch of the spiral corresponds 
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to a stable phase of the theory and that the inner branches are real instabilities char- 
acterized by a tachyonic ground state and cannot be reached by a supercooling. The 
lowest negative m branch of the spiral is tachyon free and thus could be metastable. 

The supercooling mentioned above could be attempted by considering the finite 
temperature background, namely the AdS Black hole geometry, in the presence of 
an external magnetic field. We could prepare the system in the phase corresponding 
to the trivial L = embedding and then take the T — > limit. If some of the inner 
branches of the spiral were metastable the theory could end up in the corresponding 
phase. The study of the finite temperature case is of a particular interest. Due to 
the additional scale introduced by the temperature, the theory has two dimensionless 
parameters and is described by a two dimensional phase diagram. The effect of the 
temperature is to restore the chiral symmetry and is competing with that of the 
external magnetic field. On the other side the magnetic field affects the melting of 
the mesons [50J. 

5 Acknowledgments 

V. Filev would like to thank: T. Albash, C. V. Johnson, A. Kundu and R. Rashkov 
for useful comments and discussions. This work was supported in part by the US 
Department of Energy. 



25 



A Calculating the condensate of the theory. 



Let us consider the on-shell action: 

pmax 



// F?4 /7"2 

^ V 1 + p 2 + L(p)2 Vl + ^'IP) 2 (69) 

' 

it diverges as p max — > oo. To rectify this we regularize the action by subtracting the 
action for the L = embedding: 



2n 2 T 7 N } 



pmax j 

// R 4 H 2 1 
dpp 3 Jl + - p ^ l p 2 max Vpfna X + R 4 H 2 + (70) 



l R A H 2 In ^ p2max + V R4H2 + P 



max 



4 V R 2 H 

This results to the following regularized action: 

Sreg[L, Pmax] = S — S su b (71) 

On the other-side the hamiltonian density of the theory can be written as [6]: 

H = [ d 2 9m q QQ + n , (72) 



where 7{q is the mass independent part of the hamiltonian and Q, Q are the two chiral 
fields of the hypermultiplet in M = 1 notations. Now by making the identification: 

(H) = - lim S reg (73) 

pmax *OQ 



and using equations (69)-(71) and the asymptotic of L(p) as p max — > oo: 



L{p) =m + + ... , (74) 

Pmax 



we can obtain: 



A = (#) = -2W lim S -§*L = (75) 
= lim 4TT 3 a'T 7 Nfp 3 max L\p max ) = -8n 3 a'T 7 N f c . 

Pmax ''QC 

And hence we finally get: 
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